function A=ort(n);

aux=rand(n,n);

[A d]=eig(aux'*aux);